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Abstract 

A theory of lepton decay constants based on the path-integral for¬ 
malism is given for chiral and vector mesons. Decay constants of the 
pseudoscalar and vector mesons are calculated and compared to other 
existing results. 


1 Introduction 

The decay constants f n in many cases may be directly measured in experi¬ 
ment and are important characteristics of mesons, where different theoretical 
approaches may be compared and their accuracy be estimated. The f n of 
light mesons have been studied in potential models [U [2j [3] HJ |5j 6] \ 7 \ IS, 9] in 
the QCD sum rule method [101 HU H21 [13l [lTI [151HH1 HZ] in chiral perturbation 
theory [181EEH] , as well as in lattice simulations (20] (2TJ [22] [23j [ 2H(25] [26j [27] 
and in experiment [2S1 ESI ESI El] • 

The important role of f n in theory and experiment is well illustrated by 
f n - the pion decay constant - which is the basic element and the natural 
scale of the chiral perturbation theory [32] , In the latter case f n is taken 
from experiment, and the computation of f n from the first principles is a 
serious challenge for the theory. On the lattice side a reasonable accuracy 
in computing f n was achieved recently jSnUZHEB], analytic methods include 
earlier attempts in the instanton vacuum [33] and within the Field Correlator 
Method (FCM) [35, 36 1 [37 ] , 

The present article is devoted to the systematic derivation of meson 
Green’s functions and decay constants f n for channels with arbitrary quan¬ 
tum numbers. 
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This paper is an update and extension of the earlier papers [3B], devoted 
to the heavy-light pseudoscalar and light vector mesons and [32], devoted 
to heavy-light mesons. Those papers appeared before the systematic formu¬ 
lation of FCM and in particular of the path-integral Hamiltonian based on 
FCM [50], therefore some steps in [38] |39j required a corrections. In partic¬ 
ular, the rigorous derivation of the path-integral expression for the Green’s 
function in [40] has allowed to obtain the improved expressions for decay 
constants, which we exploit in what follows. Moreover Chiral Symmetry 
Breaking (CSB) was not incorporated in [381 [39], In the present paper we 
present the consistent and general treatment of the meson Green’s functions 
and its spectral properties also for pseudoscalars accounting for CSB in the 
lowest states ( 7 t,K). The main problem which one encounters, when ad¬ 
dressing the spectral properties in QCD, is the necessity of the quantitative 
nonperturbative methods, which describe the main dynamical phenomena: 
confinement and CSB. 

The familiar (relativistic) potential model lacks the latter, while other 
models like instanton vacuum model, lacks the former QCD phenomenon. 

In what follows we are using the Field Correlator Method, which was 
introduced in [35] [36^ [37] has acquired the full form in [30] as a main tool 
to study and explain confinement. With respect to the QCD spectrum one 
derives in the FCM the effective Hamiltonian, which comprises confinement 
and relativistic effects, and contains only universal quantities: string tension 
cr, strong coupling a s and current quark masses m q . The simple local form 
of the Hamiltonian which will be called the path-integral Hamiltonian (PIH) 
occurs for objects with temporal scales larger than the gluon correlation 
length A ~ 0.2 fm, i.e. it is applicable for all QCD bound states except 
possibly toponium. 

Explicit calculations of masses and wave functions using PIH have been 
done recently for light mesons [01], heavy quarkonia [02j and heavy-light 
mesons [43], and demonstrate good agreement with experimental masses. 

In the present paper we devote a special attention to the chiral mesons and 
trat the chiral symmetry breaking (CSB) phenomenon within the formalism 
of [33], where CSB is the consequence of the confinement and the necessary 
relations can be derived from the basic parameters of QCD: string tension 
a and vacuum correlation length A, so that a fundamental quantity entering 
/tt, fk and m*-, m k is M(0) « a\ [44]. 

The paper is organized as follows: in section 2 the general path integral 
form of the meson Green’s function is presented, while in Appendices 1-5 the 
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details of derivation are given. In section 3 the obtained expression for f n is 
analyzed. In section 4 light pseudoscalar mesons are considered. The vector 
meson decay constants are studied in section 5. Section 6 contains summary 
and concluding remarks. 

2 The meson Green’s function in the path 
integral formalism 

We start with the one-body Green’s function. 

The path-integral representation for Si is [36] E23 

/s /‘OO ^ 

Si(x,y) — (rrii—D®) / ds i (Dz) xy e~ Ki $®(x, y) = (mi-D^)Gi(x, y), ( 1 ) 

J 0 


where 


1 f S * ( 
Ki = rti- Si + - j dTi ( 


dzW' 

dTi 


(x, y) = P a Pf exp ( ig / A^dzf 


x 


( 2 ) 


x exp 




(3) 


Here F^is a gluon held tensors, Pa, Pf are ordering operators, = ^( 7 ^ 7 ^— 
7 ^ 7 ^). Eqs. (PSD hold for the quark, i — 1, while for the antiquark one should 
reverse the signs of g. In explicit form one writes 


( crH crE \ 
^ crE crH J 


(4) 


The two-body q±q 2 Green’s function can be written as |20j S3j 


r 00 roc ^ 

G qi q 2 (x,y)= dsj ds 2 (Dz w ) xy (Dz^) xy {TW a (A)) A e- Kl - K \ (5) 

J 0 ^ 0 

where 

f = tr(Ti(mi — T>i)r 2 (m 2 — Do)), ( 6 ) 
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“tr” is the trace over Dirac and color indices acting on all terms. Here 
(W a {A)) is the closed Wilson loop with the spin insertions and one should 
have in mind, that color spin insertions in general do not commute, which 
should be taken into account when computing spin-dependent part of inter¬ 
action, see [l5j, in (|5]) this fact was disregarded. 


W a (A) = P a P F exp 


ig £ A^dz^ + g J cr^J F lw dr i - g J a { ^F lu /lr 2 


• (7) 


As a result of the correlator averaging 
dependent terms, one obtains 


, and neglecting the spin- 


{(W)) = Z w exp(— /Vo {r{t E ))])dt E ), ( 8 ) 

Jo 

where r{t E ) = |z x (t E ) - z 2 (t E )\, and 


V 0 {r) — V con f(r ) + Voge(i"), (9) 

rr roo 

Vconfip ) = 2 r dX duD( A, v) —> err, (r —> oo), (10) 

Jo Jo 

r oo roo 

cr = 2 du dXD{y 1 X) 1 (11) 

Jo Jo 

rr roo A ry 

V OG E= Xdx dvDf rt { = (12) 

Jo Jo Sr 


At this point it is useful to introduce as in j3D] the virtual quark (anti¬ 
quark) energies u;i(u; 2 ) instead of proper times: s t = A— , where T is the 
Euclidean time interval, T = x 4 — y 4 . As a result one obtains 


T r°° duji r°° du 2 




xy 


8vr Jo A/ 2 


n 3/2 

/0 W 2 


(D 3 0 1 ) X y(D 3 z 2 )xye- A ^^’ Zl ’ Z2 ), 


where A = K i(wi) + K 2 (u 2 ) + JV 0 (r(tE))dtE , and 


(13) 


A) (a^j) 


m , 2 + uf 


2u, 


r 1 

-T + 

Jo 


j, Ui 

dtp— 


' dzW 
dt E 


We can also introduce here the two-body 3d Hamiltonian H(coi, cu 2 , p l5 p 2 ) 
and rewrite (1131) as 
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{m\ 




xy 




( 14 ) 

where H is obtained in a standard way from the action A(a;i, cj 2 , zi, z 2 ) (we 
omit all e.m. fields except for external magnetic fields B) 


H = 


2 

E 

i=i 


p 


w ) 2 + rrij + ujf 


2l irfr 


+ Vo(r) + K,.; + AM S e 


(15) 


and Vo is given in (HU . The spin-dependent part of H, V ss and Vls are 
obtained perturbatively from cr^ u F^ w terms in (127T) . and is calculated also in 
the presence of m.f. in [45] . It is considered as a perturbative correction and 
is a relativistic generalization of the standard hyperfine interaction, 


Vss(r) = / {<r$Fi*v{x)<TpxF pX {y))d{x4, - j/ 4 )- 


Its explicit form is given in [45]. Finally, the correction —— , where 
i refers to the same quark (antiquark) yields the spin-independent self-energy 
correction AMse which was calculated earlier [46] and for zero mass quarks 
and no m.f. is 


AMse = 


3(7 

2euji 


3(7 

2huj 2 


(16) 


For the case of nonzero m.f. the resulting AMse is given in [45]. We can now 
write the total Green’s function of q L q 2 system, denoting by Y the product 
of projection operators Y = T(m 1 — Zb t )F (rn 2 — D 2 ), 


mi — Di = m 1 —ipi=mi+u 1 'y4 — ip'y, m 2 - D 2 = m 2 - cj 2 7 4 - ipj, (17) 


where p is the quark 3 momentum in the c.m. system. 
As a result one has 



y )G(x,y) = f d 3 (x — y)tr 


4Y r 


Z7T JO JO uj 2 ' 


xy 


( 18 ) 
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We have used in (fTHj) the relations 4(F) = tr(T(mi — ipi)Y{m 2 — ip 2 )), 
and neglect spin dependent terms in H\ we have taken into account, that D M 
acting on Wilson line, i.e. D M exp ('<9 J 1 A^dz^A yields exp {ig f x A^dz^d^ A. 
The c.m. projection of the Green’s function yields 


/ 



yXxle-^’^’Pi'P^ly) = X)^n(0)e _AfB(fa,1,Wa)T . 

n 


(19) 


Here M„(w 1 ,w 2 ) is the eigenvalue of H(ui, u 2 , Pl , p 2 ) in the c.m. system, 
where P = Pl + p 2 = 0; Pl = p = -p 2 . 

The integrals over dui, du 2 for T —* 00 can be performed by the stationary 
point method, namely one has 


J G(x,y)d 3 {x- y) 


T r°° du 1 
2 vr Vo 


f°° da ; 2 

, 3/2 

3 cn 2 


Y,e- Mn{uJl ’ U2)T p 2 n {0)(Y) 


where 


E 

n UJ 


e -M„( M ;°'. 4 W )T y 2 (0)<y) 

[ 0) i4 0) \/ (d 0) K (1)) (4 0) K (2)) ’ 


dM n (ui,u 2 ) 


duji 


( 0 ) 

UJ i= Uj) 


= 0 , Af”(i) = 


dM n {u ll u 2 ) 


du? 


( 0 ) 


( 20 ) 


( 21 ) 


and we have neglected the mixed terms EpE ^ or simplicity, however should 
keep them in concrete calculations: see exact result in Appendix 1 . Com¬ 
paring the results (TT9j) . (I 2 QT) with the definitions of quark decay constants 

fn 

Jr ? 


f Gr(x)d 3 x. = E / d 3x (0|jr|n)(n|j r |0)e 

J n J 

= Sr Cg) £r 

n 

where for T = 7 ^, 7^75 




d 3 P 


2 A 4 ( 2 tt ) 3 


{Mnfy) 2 ~m„T 

2 Mr, 


QnQu 


E = i, 

fc=1,2,3 y 


( 22 ) 


(23) 


and £r = 1 for T = 1 , 75 , one obtains the expression for fp (to lowest order 
in V ss ) 
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(/r) 2 = Af f r !lr:. (0)|2 . e„^\/(4 0) M;(l))(4 0) M;(2)), (24) 


, ,(°). ,( 0 ) TUT c 


This expression coincides with the previously derived in [33] [33], when 
£ n = 1/2. In what follows we show, that £ n is close to that value, but 
different for light, heavy-light and heavy-heavy mesons. 


3 Analysis of the obtained expressions 

First calculate M n , (p n from the equation 

H<p n = M n ip n , (25) 

treating Vls, Vss and Vse in (9) as perturbation, H = + Vls + Vss + 

A Mse = + AM n . One simplify the procedure introducing the relative 

coordinate in the c.m. system, rj = iq — r 2 , p = so that without magnetic 
field 

r\2 _ TfP -4- t,p 

= ^ + + V o(r) + A M se , (26) 

2u 2oji 

H = H 0 + V SS + V LS . 

(Note, that this is not nonrelativistic expansion!) and finding stationary 
values of M^\ H 0 ip n = M^ip n with respect to aq,aq = ) from the 

equation 



duji 


,( 0 ) 


0 . 


(27) 


This is the basic approach in the string Hamiltonian formalism m and 
it was checked that the accuracy of the replacement aq = for lowest 
states is around 5% m ■ However the values <^ n (0) are more sensitive to the 
replacement (23), and one should use original Hamiltonian (15) to calculate 
(p n ( 0) [38] 09], see Table 4 of j38] for comparisons. 

It is essential, that we are using Vss + Vls as perturbation terms to 
compute the final hadron masses and hence M n in different parts of ([24]) is 
inserted as computed from H 0 ([26]) . not containing Vss + Vls- 

ii) As it was argued above, the factor Ip can be computed in terms of 
momenta of quark and antiquark, or in the c.m.system in terms of relative 
momentum p, with the result. 
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Y v = mim 2 + oj itu 2 + -p 2 , (28) 

o 

Y S = —17111712 + UiQ 2 \ + P 2 , (29) 

p 2 

Y Ai = -17171712 + UllUl2 + —, (30) 

Ya 4 = m 1 m 2 + LO 1 UI 2 ~ p 2 , (31) 

Y P = (mim 2 ) + u>iu) 2 - p 2 , (32) 

Here we used notations: 

Y v = \ tr [( m \ - Di)li{m 2 - £> 2 )7i], (33) 

i 6 

Y Ai = - Di)^a 5 (m 2 - £> 2 ) 7 * 75 ], (34) 

d i 

Ya 4 = -[trim! - 1 ) 1 ) 7475(^2 - ^ 2 ) 7475 ]• (35) 


In case of the pseudoscalar channel in (j3TT) . (j32j) in the chiral limit 
mi,m 2 —> 0 there appear additional mass terms due to CSB, which are 
computed through held correlators and are given in [Hj. As it is shown in 
Appendix 2, the proper account of CSB leads to the fact, that Eq. (1241) for 
fp retains its form for it, K mesons, but the expression for chiral mesons e.g. 
Ya 4 should be replaced by more general one, 

Yj 4 ^ = (mi + Mi (0)) (m 2 + M 2 (0)) + U 1 U 2 — p" (36) 

In the chiral limit mi = m 2 = 0 it was found in |33| that Mi(0) = M 2 (0) = 
0.15 GeV and it was computed through the held correlators, M(0) = aA, 
where A is the vacuum correlation length, A ~ 1 GeV^ 1 [0Hj. 

In the nonrelativistic limit, m* y r a, one can easily find that uj, fa mi, 

while (p 2 ) ~ 0(a) , and therefore one has 

(Yv)nr ~ 2mim 2 + 0(a), ( Y s )nr ~ 0(a), (37) 

(Ya^nr ~ 0 (a); (Y Aa )nr = 2mim 2 + 0 (a), (38) 


{Yp)nr = 2mim 2 + 0(cr). (39) 

Therefore in the nonrelativistic limit m\ yfo, m 2 y/a, for f n for V and 
P channels one obtains 


fn(NR) 


4 Nc 
M n 


kn(0)P 


(40) 


as was found earlier [1], 

As a final step one needs to compute the radiative corrections to f n , 
which come from the short-distance (large momentum) perturbative gluon 
contributions. Neglecting interference terms they can be written as in [ 23l 

El!, 

(W a ) = ( W OGE )(W nonpert ) (41) 

and 

, TTr , ^ / 4 /• r dz±dz' A a s (z — z')\ , 

(W oaE ) = exp (--]] ) , (42) 

where Z m is a regularization factor. After separating the Coulomb interaction 
in H in this way, one gets the correction to {W a ), and /f can be written as 

,/r “^ ^r/n ■Cr = 1 + c r a s + 0(a 2 s ). (43) 


but this correction is small and will be neglected below. 

Another important contribution from perturbative gluon exchanges (GE) 
is the account of the running coupling constant in (40) which is especially 
important for ip n ( 0) in the S'-wave channels. Introducing the asymptotic 
freedom factor Paf (we follow notations from (38] ) 


Paf 


y^ F) (Q) 

^n(O) 


(44) 


one should multiply with this factor and finally gets 

fr = fr irPAF- 


(45) 


We conclude this section with the discussion of input parameters in the 
approach described above. The set of parameters includes m.j, a s and a in 
the first approximation and mj, a s (r), a, when asymptotic freedom is taken 
into account. 
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Here m* are pole masses which are connected to the Lagrangian (cur¬ 
rent) masses in MS scheme as (see [28] for review and [39] for additional 
references). 



4 Light pseudoscalar mesons and current cor¬ 
relators in PS channels 


The formalism of the previous section is of general character and be applied 
to any channel T. However to save space we shall consider below only PS 
and vector mesons. Pseudosclalar mesons appear both in A 4 and P channels. 
Their connection to the A channel is given by the chiral anomaly term: 


(0|j£ A) |g,n) = i/p(n)^. (47) 

Exploiting this definition one obtains the same expression as before for 
fp(n), considering the c.m. system P = 0, namely 


H(n) 


V(df r °' ) )ly»( 0) l 2 


(48) 


One should have in mind however that our formalism above in this paper 
till now did not take into account Chiral Symmetry Breaking (CSB) and 
therefore cannot be applied to the Nambu-Goldstone mesons 7 r, K, 77 . For the 
latter one should use the technic suggested and exploited in [ 33 ], where 
was computed through the masses m n , <^ n ( 0 ) as in (t24|) but in addition there 
appears an effective mass parameter M(0), see Appendix 2 . The resulting 
equation for (Ta 4 ) Eq. pfil) can be written in the limit m, —> 0 as 

(Ya 4 ) = M 2 ( 0 ) + (^/p 2 + mf)( v /p 2 + m\) - (p 2 ) M 2 ( 0 ). (49) 

I 11 the chiral limit, rn\ = m 2 = 0 , and taking into account that Oi = 
uj 2 = and (p 2 ) = oj 2 , one has 

fUn) - (50) 

(jJ i\// n ^ n 
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For <p n (0) and M n one takes neglecting hyperfine interaction the same 
values, as for p-rneson, i.e. M n = M(n = 0) = 0.65 GeV, |<p n (0 )| 2 = 
0A09 £ eV3 , u = 0.352 GeV. 

Taking now = 2.45 from A(12) and M(0) = 0.15 GeV one obtains 
/ 2 = 0.01782 GeV 2 , f n = 133 MeV. This should be compared with the 
experimental value, which in the normalization of Eq. (I5U1) is equal to f* x = 
\pl ■ 0.093 GeV = 0.131 GeV. One see that agreement within 2 %. 

One should stress, that in absence of CSB, when M(0) = 0, and (Va 4 ) —> 
m\rri 2 —> 0, also f n vanishes, implying that f n plays the role of the CSB order 
parameter (together with (qq), which is also proportional to M( 0 )). 

We now turn to the case of K meson. Doing calculations in the same way 
as for pion above, and taking m s = 0.15 GeV, a = 0.18 GeV 2 , one has for 
K- meson; 


l u u {K) = 0.36GeV, ou s (K) = 0.39GeV, M^ 0) = 0.84 GeV. 

The latter number is obtained without Coulomb and hyperfine interac¬ 
tion, which shift M k by A mcoui = 0.05 GeV and A muf = 0.3 GeV, resulting 
in rriK = — A m cou i — A rriHf — 0.49 GeV. 

Using Appendix 1, Eq. (A2), one obtains ^ = 2.29 so that f\ is 



2.29 • N c (Y K )<p 2 K (0) 

uj u uJsMk } 


(51) 


where (Y K ) = (M(0) + m u )(M( 0) + m s ) + oj u ujs — (p 2 ) = 0.06 GeV 2 and as 
a result Jk = 0.165 GeV. 

To compare Jk and f n we write down for both mesons without hyperfine 
interaction and using (49) 

% = 1-6, (52) 

fn 

tj- = 1.24 

Jn 

The result f)52|) . = 1.24 is in agreement with the experimental values 

J 7T 

m 


= 130.7 ± 0.1 ± 0.36 MeV, 


(53) 
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Table 1: Decay constants of chiral mesons and its excitations 


7r (nS) 

K(nS) 

n 

1 

2 

3 

1 

2 

3 

/™(GeV) 

0.138 

0.069 

0.048 

0.165 

0.104 

0.085 

fwn/ fir 1 

1 

0.5 

0.35 

1 

0.63 

0.515 


/^ p) = 159.8 ± 1.4 ± 0.44 MeV, 

which yields 


/ (exp) 

I K+ 


/, 


(exp) 


= 1.22 ± 0 . 02 , 


(54) 


while lattice calculation [20] yield for this ratio 1.195 ± 0.006. 

We now turn to the radial excitations of the chiral mesons. In this case 
of high excitations there appear decay channels, which play a role of inter¬ 
mediate channels in the meson Green’s function. Therefore neglecting these 
channels, we shall make only rough upper limit estimates of decay constants. 
To this end we exploit the fact (see Appendix 1) that £ n does not depend 


on n, and <^(0) can be estimated to the lowest order as 


Xn(x) 


M n « 


X)n(O) 

2oj n ( with the accuracy of (10± 15)%). As a result one obtains 

0.105 


Un( GeV) ^ 


M n { GeV)’ 


n > 1 


and 


(55) 


where M n are the mass values before the chiral shift |49], and as a result one 
has the values listed in the Table 1. 

One can see, that f nn and fe n are slowly decreasing for growing n, im¬ 
plying a substantial leptonic decay contribution to the list of decay modes. 


5 Decay constants of vector mesons 

We now turn to vector mesons, where we take the same value for M n and 
£n (H , not affected by the hf splitting, which we take afterwards as a 
perturbation. Hence we take fe 1 = 2.29 for K* meson and = 2.45 for p 
and oj mesons, while for the 0 meson one obtains £0 1 = 2.095, and cu s = 0.424 
GeV, A M S e = —0.282 GeV and = 1.040 GeV. (M^(exp) = 1.020 MeV). 
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Table 2: Decay constants of vector mesons 


Vi 

P 

UJ 

0 

/■!.„ ( GeV) 

0.254 

0.0846 

0.096 

/ 17 (exp) (GeV) 

0.255 

0.0756 

0.107 

T ee ( ex p) (keV) 

7.04 ± 0.06 

0.60 ±0.02 

1.26 ±0.02 


Since vector mesons are connected to the electromagnetic current, the 
corresponding decay constants contain the effective charge e 2 (i) = -V, | 

for % — p, uj, 4> [] and vector decay constants have the form 



e 2 q m^imyy) 

wi 0) o4 0) M n (i)€„(i) 


P,u,(p, 


(56) 


while the dielectron width is connected to f^\i) as [28] 


T(V -)• e + e~) 


Ana 2 

3M v (i) 


C/f’M) 2 ) 1 



(57) 


Keeping the value of the last factor in (1371) to be equal 0.32(a s = 0.4) 
for p, uj and cq = 0.3 for <p, as a result one obtains the values of / v (exp), 
given in the Table 2. Lattice data [22] yield f p = 239(18) in a satisfactory 
agreement with our result in Table 2. 


6 Summary and conclusions 

We have presented the theory of lepton decay constants for light mesons 
based on the path integral formalism. It essentially exploits the path integral 
Hamiltonian (PIH) depending on virtual q, q energies UJ\,U 2 and the final 
form obtains after the stationary point analysis with respect to U \, u- 2 - The 
same approach has given a large number of observables (masses and wave 
functions, Regge trajectories etc.) in good agreement with experiment both 
for light and heavy quarks munis]. The lepton decay constants for heavy- 
light meson have been calculated in the same method in an approximate form 
in [38], 39] also in good agreement with lattice and experimental data. 

In this paper we specifically considered light mesons and paid a special 
attention to the chiral mesons and their radial excitations. We have also 
exploited an improved form of the path integral from [40], which allows to 
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obtain much better accuracy. Another important ingredient of the present 
paper is a new treatment of chiral mesons, which exploits the fundamental 
quantity - the scalar chiral mass parameter M( 0) ~ a\ (the corresponding 
chiral correlation length is 1/crA), which as shown before in [33] and here in 
Appendix 2, enters additively with the current quark mass m q and disappears 
at large m q . 

The resulting values of f n and fx are in good agreement with experimen¬ 
tal data, however the only parameters of our theory are rn q , a s and string 
tension a. We have also calculated decay constants of vector mesons p,u,(j) 
and found a satisfactory agreement with experiment. 

The author is grateful to A.M.Badalian for useful suggestion and dis¬ 
cussions. Financial support of the RFBR grant 1402-00395 is gratefully ac¬ 
knowledged. 


Appendix 1 


The correction coefficient 


As it was shown in [40] (see also appendix of [30]). in (103]) is defined 
as follows: 



(Al.l) 


with 



(A1.2) 


and 


a = - 


1 d 2 M n = idfAT = d 2 M n 

2 duo 2 ’ 2 diij\ ’ ^ diij\diijT 


(A1.3) 


M n is defined as 


i=l,2 



_I uji ~r u. 

cu =- 

U1UJ2 


UJi + U 2 


(A1.4) 
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and e n is the eigenvalue of the equation 


P1 

2 Cj 


+ V(r) (p n (r) = e n (p n (r), 


(A1.5) 


where V(r) includes confinement V c (r) and gluon exchange interaction V g (r), 
but not hyperhne interaction, which is taken into account as a first order 
correction to the total mass, and hence is not to be present in (1A1.3D . Hence 
finally M n , entering in f 2 , is the hyperhne averaged eigenvalue M n with the 
selfenergy term A$e taken into account 


v„ = y (lJ,(o y,ot m ‘ + e „((* <0) )- 1 ) + a se 


2 aU 


(Al. 6 ) 


For a, (3 ,7 one obtains 


2a, 2/3 — —j H- 3 e' n H— j£ n , i — 1,2. 

ujf uf ujf 
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(A1.7) 

(A1.8) 


l-Tp'rp p — p — 

neie t n — o~ 1, t n — 


d 2 e 


dti- 1 ’ n 0 (d)— 1 ) 2 ' 

For mi = m 2 = m and hence oji = u >2 = to one has a — (3 and is 


n _ 2 7 2 _ £ n( £ n + If) + + If)" 

“ "n ™ ^ — 


4cu 7 


In the nonrclativistic limit e n <^im,uj = m, one has 

n — 1 t — ~ 

71 Am 2,U 2 


(A1.9) 


(A1.10) 


Consider now light quarks and put m\ = m 2 = 0, and first neglect the 
OGE interaction. Then 


£„ = 2 -‘/3( i ,-.)l/3 (7 2/3 o(n ) ! 

with ao = 2.338 and other a(n) given in Table 2 of [47] , 
Using dsrzD, CM) , one obtains 


£ n (a s = 0 ,m; = 0 ) = 


= 0.408; C 1 = 2.45. 


a/54 


(Al.ll) 


(A 1 . 12 ) 
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The resulting / 2 is 


f 2 - 
Jr.n — 


2.45iV c (lr)^ 2 (0) 

oJoM n 


where <^ 2 (0) = ^, (Y v ) = |wq> M) = 4w 0 + A SE = 2w 0 . 
The inclusion of OGE interaction yields ESI S3 

4 S) = = 2- 1 /3( l5 - 1 ) 1 /3^/3 o(A 


a(n 


(A1.13) 


(A1.14) 


where A = ^ ^ and a(0, L, n) = a(L,n). 

The values of a(A, 0, 0), ff(A, 0, 0) are given in the Table 4 of p3j, and 
one has an estimate of |^f = 0.2 for A < 0.9. 

The most important change due to OGE is in </? 2 (0) which is now for 


L = 0 


¥>n(°,«-) = + ^a{r 2 )n) 


<7UJ_ Xa(0) 
4tt Xq(0) 


(A1.15) 


The values of </? 2 (0, a s ) are given in the Table 4 of [38], together with 
xa(Q) 2 

Xo(0) 


Appendix 2 

Chiral correction length in the confining vacuum 

It was shown in [44] . that nonzero field correlator (FF) = (tr(F(u)<j)F[v)<p ), 
generating the kernel J(x, y) ~ f du f dv(FF), leads to the appearance in the 
quark Green’s function S(x,y ) the nonpertnrbative nonlocal mass operator 
M(x,y), M(x, y) JS, satisfying equation 

0 + m q )S(x, y) + J M(x,y)S(z,y)d A z = h 4 (x - y). (A2.1) 

This general property can be further analyzed separating in A4 scalar-isoscalar 
part M. s and pseudoscalar-isovector pieces U which can be conveniently writ¬ 
ten as 
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M —> MsU, U = exp(?750) (A2.2) 

As a consequence the scalar nonlocal mass A 4s(x,y) enters into the quark 
Greens function S' as a scalar piece together with the current quark mass m q , 

iS(x,y) = 0+ m q +M s )~y. (A2.3) 

Note, that at large m q , the magnitude of M s is fast decreasing 

M s ~ 0(l/m q ), m q — > oo. (A2.4) 

In the current correlator (x|rS' (? rS'q||/), Ms enters in the local vertex form 
Ms(x, x ) = M s ( 0), and in the framework of the chiral approach |44| based on 
(IA2.1() and (1A2.2|) . one derives the relation for in the chiral limit, m q —>- 0, 

ft = N C M\ 0) £ (A2.5) 

n =0 m n 

where n refers to the radial excited qq states with mass M n , and the sum 
over n is cut off by the factor exp(— M n X),X = 0(1 GeV -1 ) is the vacuum 
correlation length, calculated via gluelump masses [38], M(0) was calculated 
in the Appendix 4 of the third reference [43], 

Ml 0) = « 0.15 GeV. (A2.6) 

\Ar 
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